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Cn i Abstract. We show that two qubits initially in completely classical state can create 

quantum discord through a local generalized amplitude damping channel, but high 
rS^ ' temperature will impede the creating of quantum discord. 
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1. Introduction 

^ ■ Quantum entanglement is one of the most striking features of quantum physics. It 

t^^ ■ remarkably reveals nonlocality and may lead to powerful applications in quantum 

ly-^ ■ information and quantum computation [H |2]. But entanglement is incredibly fragile 

f~^ ■ under environmental noises, and any local operations can not increase it. These are 

serious frustrations for using entanglement. 

Quantum discord [Sj H] is another quantum correlation other than entanglement, 

and also can be used for quantum computation [SIEIE]- It has been shown that quantum 
^ . discord is more robust than entanglement [8]. Furthermore, recent studies [9l [TOl [11] 

indicate that local environmental noises can create quantum discord. These are exciting 

results compared to the case of entanglement. 

In [10], F. Ciccarello and V. Giovannetti showed that a zero discord state of two 

qubits 



O 



S 



r,n. = ^(|+)(+|®|0)(0| + |-)(-|®|l)(l|) (1) 

can create quantum discord under an amplitude damping channel, where 

l±) = ^(|0)±|l)). (2) 

In this paper we will investigate the more general case that the initial state is 

P^n^ = ^(|^o(A)) (^o(A) | ® |0)(0| + |^i (A)) (^i(A) | ® |1)(1|), (3) 

where 

\M>^)) = V\\o) + VT^^\i), \M>^)) = VT^\\o)-v^\i), AG [0,1], (4) 
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and the noise is the generahzed amphtude damping channeL We know generahzed 
amphtude damping describes temperature effects, so we are especially interested in the 
temperature effect on the quantum discord. This paper is organized as follows. In 
Sec. 2, we briefly recall some basics about quantum discord and generalized amplitude 
damping as preparations. In Sec. 3, we investigate the quantum discord of two qubits 
whose intial state as in Eq.(3) and undergoes a local generalized amplitude damping. 
Sec. 4 is a brief summary. 

2. Quantum discord and generalized amplitude damping 

Suppose two quantum systems A and B are described by the Hilbert spaces H"^ and H^ , 
respectively. The composite system AB is then described by the Hilbert space H^® 
H^. For a state p on H"^® H^ , the quantum discord (with respect to A) of p is defined 

as [3] 

DM = Sip^) - S{p) + inf j:P^S{pf), (5) 

where p^ = trsp, Hj = \i){i\, inf takes all projective measurements {IIjjj on system A, 
Pi = trB{i\p\i), pf = {i\p\i) IVi) S{-) is the Von Neumann entropy. It can be proven that 

m 

Da{p) =0^^p = Ep.K)(^I ® pf, (6) 

i 

where {\i)}i is an orthonormal basis for H^, {pf}i are density operators on H^, pi > 0, 

EiPi = 1- 

Eq.(5) is hard to optimize, even for two qubits case till now only few states are 
found possessing analytical expressions [121 113]. B. Dakic, V. Vedral and C. Brukner 
proposed a geometric measure of quantum discord as [H] 

Di{p) = \ni{tr[{p-af]:DM) = ^}- (7) 

Evidently, 

Da{p) = ^^ Diip) = 0. (8) 

As an elegant result, D^[p) allows analytical expressions for all two qubits states [H]. 
More specifically, for a two qubits state 

13 3 3 

p = -{I ® I + Y^ XiCJi (g)i + J2yj^ ® ^j + Yl ^ii^i ® ^i)' (9) 

4=1 j=l i,j=l 

where, cxi = a^, o"2 = o"?/! ^"3 = a^ are Pauh matrices, {xjjf^^, {yj}^^^, {Tij}\j^i are all 
real numbers with 

Xi = tr[pai ® I],yj = tr[pl (g) cXj],Tij = tr[pai (g) ctj], (10) 

then [H] 

^^(P) = l(E^'+E^?-Vax). (11) 

^ i=l ij=l 
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Where Amax is the largest eigenvalue of the matrix xx* + TT^, x = (xi, X2, X3)*, T is the 
matrix (Tj^). 

Generalized amplitude damping describes the effect of dissipation to an environment 
at finite temperature ([1], 8.3.5). Suppose two qubits systems A and B, given the bases 
{|0), |1)}, {|0), |1)} for H^ and H^, the operation elements for generalized amplitude 
damping are 



Where p G [0,1], 7 G [0,1]. p = p{T) and 7 = 7(T) are functions of temperature 
T. 7(T) describes the transition probability between |0) and |1) at temperature T. 
p{T) describes temperature effects, p(0) = 1 is the case of amplitude damping. We can 
assume that (see for example [I5]) p{T) is a monotonically decreasing function when T 
varies from to +00, and 

limT^opiT) = l,limT^+ooP{T) = 0.5. (13) 

From Eq.(12) we can see that p varies from to 0.5 is symmetric to that p varies from 
1 to 0.5 in the sense exchanging states |0) and |1), so we only need to consider the case 
p varies from 1 to 0.5. 

A state p after this generalized amplitude damping will be a state 

E{p) = EopEl + E.pEl + E^pEl + E^pEl (14) 

3. Creating quantum discord through local generalized amplitude damping 

Now suppose two qubits initially is in the state as in Eq.(3), through the generalized 
amplitude damping channel in Eq.(12), the output state p = E{pini) can be calculated 
directly by Eqs.(3,12,14), will be 

p=-{I®I + X3(T3 (S)I + Tisai ® (73 + Tss^a ® ag), (15) 

where 



X3 = (2p-1)7, Ti3 = 2v/A(l-A)(l-7), T33 = (1 - 7)(2A - 1). (16) 

We first ask, for what initial states and what generalized amplitude damping channel, 
the quantum discord of the output states, Da{p), are vanishing? This is the Proposition 
1 below. 

Proposition 1. The quantum discord Da{p) of the state p in Eq.(15) is zero if one 
of the conditions below is satisfied 

(i). 7 = 0; 

(ii). 7 = 1; 

(iii). A = 0; 

(iv). A = 1; 
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(v). j9 = 0.5. 

Proof. From Eq.(8), we only need to prove D^{p) = 0. Using Eq.(ll) to Eq.(15), 
after some calculations, we get 

D^ip) = \{i4 + T^3 + Tl) - v/[(a:3 + Ti3)2 + T|3][(x3-Ti3)2 + T|3]}. (17) 

From Eq.(16), it is easy to see that 

(vi). if 7 = 0, then 0:3 = 0; 

(vii). if 7 = 1, then T13 = T33 = 0; 

(viii). if A = 0, then T13 = 0; 

(ix). if A = 1, then T13 = 0; 

(x). ii p = 0.5, then X3 = 0. 

Each of the cases above can result in D'^{p) = 0. We then complete this proof. 

We next investigate the quantum discord Da{p) for the state p in Eq.(15) for more 
general cases, that is 7 G (0, 1), A G (0, 1), p G (0.5, 1]. 

Any orthonormal basis of H^ can be expressed as an unitary matrix U multiplied 
by the orthonormal basis {|0), |1)}. Since U can be written as 

\h + ia t-ic ) ^ ' 

where t, a, 6, c are all real numbers and satisfy 

t^ + a^ + 6^ + c^ = 1, (19) 

then any orthonormal basis of H^ can be expressed as 

|v9o) = (t + ^c)|0) + (6 + m)|l), 

|v?i) = (-6 + m)|0) + (t-ic)|l). (20) 

Notice that 

/|0) = |0), a,|0) = |l), ay\Q) = t\l), a,|0) = |0), 
/|1) = |1), a,|l) = |0), ay\l) = -i\Q), a,\l) = -\l). (21) 

Denote poPo = (v^oIpIv'o) and pipi = {(pi\p\(pi), then from Eqs.(15,20,21), we have 

PoPo = ^-[l + it^ + c^-a^- b')xs] + ^[2(t6 + ac)Tis + {^ + c^ - a^ - b^)Tss], (22) 

Po = tr{Lpo\p\(po) = -(1 + 3:3) - (a^ + b^)x3, (23) 

PiPi = ^-[l + {a^ + b'-t^- c')xs] + ^[-2{tb + ac)T^s + (a^ + b^ - t^ - c'')Ts3], (24) 

pi = tr{ipi\p\ipi) = -(1 - 3:3) + (a^ + b^)xs. (25) 

Let 

a = y/rcos9i, b=^/rsm9l, c = \^l — r sin 62, t = \^l — r 00562, (26) 

then 

a'^ + b'^ = r, tb + ac= Jr{l -r) sin 6*3, (27) 
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where 

re [0,1], ^iG[0,27r), ^2G[0,27r), ^3 = ^1 + ^2- (28) 
Eqs. (22-25) then read 

popo = ^[1 + (1 - 2r)x3] + ^[2y^r(l-r) sin ^3^13 + (1 - 2r)T33], (29) 

Po = ^[l + (l-2r)x3], (30) 

PiPi = ^[1 - (1 - 2r)x3] + ^[-2y^r(l-r) sin ^3^13 - (1 - 2r)T33], (31) 

Pi = \[l-{l-2r)x,]. (32) 
Let 



a = 1 - 2r, 13 = 2^r(l - r) sin ^3, (33) 

then Eqs. (29-32) become 

PoPo = ifl + '^^3] + -f[PTi^ + "^33], (34) 

Po = 2[1 + '^^3], (35) 

Pipi = t[1 - "^3] + -^[-PTi^ - "^33], (36) 

Pi = -[l-ax3], (37) 

with 

a^ + (3'^<l. (38) 

From Eqs. (34-37), it is easy to find that po has two eigenvalues 

l + a(a:3 + T33) + /3Ti3 ^ _ l + a(a:3 + T33) + /3Ti3 

2(l + aa:3) ' 2(1 + aa;3) ' ^ ^ 

and pi has two eigenvalues 

1 - a{x-i + T33) - /3Ti3 ^ _ 1 - a{x-i + T33) - /3Ti3 
2(1 — aa;3) ' 2(1 — ax^) 

Thus 

PoSipo) +PiSipi) = F{a,l3) 

1/1 , . .l + a(a:3 + T33) + /3Ti3 1 1 - a(a;3 + T33) - /3Ti3 

= -l + aa;3/i -— — ^ ) + -{1 - ax3)h{ — ^ , 41 

2 2(1 + ax^) 2 2(1 — ax^) 

where h{-) is the binary entropy 

h{x) = —xlog2X — {1 — x) log2(l — x), (42) 

h{x) is defined on a; G [0, 1], h{x) is a concave function, and h{x) = h{l — x). 

We say F{a, /3) is a concave function in (a, /3), this can be seen from the concavity 
of h{x) and the facts below: 
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(xi). If f{x) and g{x) are convex functions, then so is f{x) + g{x) ([TB], 3.2.1); 

(xii). Convexity is invariant under afiine maps: that is, if f(x) is convex with x G -R", 
then so is giy) = f{Ay + Xo) with y G R"^, Xq G -R", A annxm real matrix ([I6j, 3.2.2); 

(xiii). If f(x) is convex, then so is g{x,t) = tf{x/t) with t > ([IS], 3.2.6). 

It follows that the minimum of F{a, P) in Eq.(41) over the domain a^ + /3^ < 1 can 
be achieved on the unit circle a^ + /3^ = 1. Further, from Eq.(41), it is evident that 

F{a,f3) = F{-a,~l3), 



hence, 



where 



min[po5'(po) +PiS{pi) : all {Hi} 



i=i\ 



mm{F{9) : < < vr}. 



(43) 
(44) 



F{e) 



l.T , 1 + (x3 + T33)cos6' + Ti3sin6' 



4^ 



x^cos 9)h{ 



1 - (x3 + T33) cos 6 - Ti3 sin 9 , 
2(l-X3Cose) 



Finally, we get the expression of D{p) as 

D{p) = S{p^) - S{p) + mm{F{e) : < ^ < tt}, 
where S^p"^) and S{p) can be calculated directly by Eq.(15) as 



S{p^) = h{ 



S{P) 



1 



1 + X3, 

-M- 



I,l + J{x3 + Tssy + Tl,, 



2 ^ 



(X3 - T33)2 + Ti 



13, 



(45) 
(46) 

(47) 
(48) 



2 ' 2 

From these expressions, we have Proposition 2 below. 

Proposition 2. Quantum discord Da{p) in Eq.(46) for the state p in Eq.(15) has 
the same value for (p. A, 7) and (p, 1 — A, 7), that is, Da{p) is symmetric about A = 0.5. 

Proof. From Eq.(16), if (p. A, 7) generates (X3, T13, T33), then (p, 1 — A, 7) generates 
(0:3, Ti3, -T33). (x3, Ti3, T33) and (0:3, T13, -T33) generate the same value for 5'(p^)-S'(p), 
see Eqs. (47-48). Also, in F{9) of Eq.(45), by use of h{x) = h{l — x), we can get 

F{e, X3, Ti3, T33) = F{n - 9, X3, T13, -T33). (49) 

Thus, we can readily attain Proposition 2 and end this proof. 

p=l 

^ p=0.8 
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p=0.6 



p=0.55 




O.OOD: 
O.OOOD t- 



Figure 1. Discord versus A and 7 when p = l,p = 0.8, p = 0.6, p = 0.55. 



Fig.l depicts Da{p) as a function of (A, 7) 
at p = 1, p = 0.8, p = 0.6, p = 0.55, 
respectively. From Fig.l we see that, when 
A is close to A = 0.5 and p close to 1, 
Da{p) is relatively large. That is to say, more 
superposition of the states |0) and |1) in Eq.(3) 
can boost the creating of quantum discord, 
while high temperature, conversely, can impede 
the creating. Fig. 2 shows Da{p) as a function 
of (p, 7) at A = 0.5. 

4. Summary 



A=0.5 




Figure 2. Discord versus p and 7 when 
A = 0.5. 



We investigated the quantum discord of two qubits which initially in completely classical 
state then experienced a local generalized amplitude damping channel. We showed that 
a completely classical state can create quantum discord through a local generalized 
amplitude damping, while high temperature will impede the creating of quantum 
discord. 

Acknowledgements 

This work was supported by National Natural Science Foundation of China (Grant Nos. 
10775101). The author thanks Qing Hou for helpful discussions. 



References 

[1] M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information (Cambridge 

University Press, Cambridge, UK, 2000) 
[2] R. Horodecki et al. Rev. Mod. Phys. 81 (2009) 865 
[3] H. OUivier and W. H. Zurek, Phys. Rev. Lett. 88 (2001) 017901 



Creating quantum discord through local generalized amplitude damping 8 



[4] 

[5] 
[6] 

[7] 

[8] 

[9] 

[10 

[11 
[12 

[13 

[14 
[15 
[16 



L. Henderson and V. Vedral, J. Phys. A 34 (2001) 6899 

A. Datta, A. Shaji and C. M. Caves, Phys. Rev. Lett. 100 (2008) 050502 

B. P. Lanyon, M. Barbieri, M. P. Almeida and A. G. White, Phys. Rev. Lett. 101 (2008) 200501 
A. Datta and S. Gharibian, Phys. Rev. A 79 (2009) 042325 

T. Werlang, S. Souza, F. F. Fanchini and C. J. Villas Boas, Phys. Rev. A 80 (2009) 024103 

Z. Merali, Nature 474 (2011) 24 

F. Ciccarello and V. Giovannetti, (2011) a7-Xw:1105.5551 [quant-ph] 

A. Streltsov, H. Kampermann and D. Bruss, (2011) arXit):1106.2028 [quant-ph] 
S. Lug, Phys. Rev. A 77 (2008) 042303 

M. Ah, A. R. P. Ran and G. Alber, Phys. Rev. A 81 (2010) 042105 

B. Dakic, V. Vedral and C. Brukner, Phys. Rev. Lett. 105 (2010) 190502 
R. Srikanth and S. Banerjee, Phys. Rev. A 77 (2008) 012318 

S. Boyd and L. Vandenberghe, Convex Optimization (Cambridge University Press, Cambridge, 
UK, 2004) 



